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Classical field theory is considered as a theory of unparametrized surfaces embedded in a con¬ 
figuration space, which accommodates, in a symmetric way, spacetime positions and field values. 
Dynamics is defined by a (Hamiltonian) constraint between multivector-valued generalized mo¬ 
menta, and points in the configuration space. Starting from a variational principle, we derive local 
equations of motion, that is, differential equations that determine classical surfaces and momenta. 
A local Hamilton-Jacobi equation applicable in the field theory then follows readily. The general 
method is illustrated with three examples: non-relativistic Hamiltonian mechanics, De Bonder-Weyl 
scalar field theory, and string theory. 


I. INTRODUCTION 


In non-relativistic mechanics, the trajectory of a particle is a function x{t), which describes 
how the position of the particle changes with time. In relativistic mechanics, space and time are 
treated equally, and the particle’s trajectory is regarded as a sequence of spacetime points {t, x). 

In field theory, the field configuration is usually viewed as a function 4>{x) that describes how 
the field varies from point to point. However, general relativity suggests [T] that the spacetime 
is a dynamical entity, which should be put with fields on the same footing. Mathematically, 
instead of a function (j){x) one should therefore consider the respective graph, i.e., the collection 
of points (a;, (jf). 

In this article, we study the mathematical formalism proposed in pQ Ch. 3] that treats time, 
space, and fields equally. All these entities are collectively called partial observables, and together 
they form a finite-dimensional configuration space. Classical field theory predicts that certain 
correlations between partial observables can be realized in nature. These are then called physical 
motions, and have the form of surfaces embedded in the configuration space. 

Our dynamical description utilizes multivector-valued momentum variable, which can be 
thought of as conjugated to the motion’s tangent planes; thus generalizing the canonical mo¬ 
mentum conjugated to the velocity vector in classical mechanics. Individual theory is specified 
by a choice of the Hamiltonian H, which is a function of a configuration space point q and mo¬ 
mentum P. This Hamiltonian enters into a variational principle (Section|H| via the Hamiltonian 
constraint H{q,P) = 0. 

The aim of this article is to establish, in the first place, equations of motion that follow from 
the variational principle. This is done in Section HI Eqs. (13). These equations generalize the 


Hamilton’s canonical equations of motion of classical mechanics. From Eqs. (13) we derive the 


local Hamilton-Jacobi equation (20), which generalizes to the field theory the respective concept 
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from classical mechanics (in this context, see also Refs. [ 2 ] and 0 ). It should be stressed that 
both, Eqs. (13) and (20), contain only partial, not variational, derivatives. 

Three examples are provided in Section [V] to demonstrate universality of the present formal¬ 
ism. The first example (VA) shows how non-relativistic mechanics is deduced when we assume 
that the motions are one dimensional curves, and choose the Hamiltonian H appropriately. Eqs. 
(13) then reduce to the Hamilton’s canonical equations, accompanied by the law of energy con¬ 
servation. The Hamilton-Jacobi equation of classical mechanics is also recovered. 

The second example (V B) discusses the theory of real one-component scalar field defined on 
a Euclidean spacetime of any dimension. It is shown that Eqs. (13) produce the De Donder- 
Weyl equations HE], and, at the same time, they incorporate the continuity equation for the 
energy-momentum tensor. Hamilton-Jacobi equation reproduces the one invented by Weyl [Sj. 

In the last example (VC) we treat relativistic particle, string, or higher-dimensional membrane, 
depending on the dimensionality of the motions. The configuration space is identified with the 
target space of the string theory, motions are the worldsheets, and the corresponding Hamiltonian 
is essentially the simplest and most symmetric function of the momentum variable. The equations 
of motion have simple geometric meaning, namely, they ensure that the mean curvature of the 
physical motion vanishes. In fact, this is exactly the condition that defines minimal surfaces [ 8 |. 

One more remark is in order before we start. All manipulations are performed in the mathe¬ 
matical formalism of geometric (or Clifford) algebra and calculus developed by D. Hestenes [Sj. 
It is a coordinate-free language that is more universal than the calculus of differential forms, 
nevertheless, it is yet not well-recognized by a broad audience. Reader unfamiliar with geometric 
algebra or calculus is recommended to first read Appendix]^ where we introduce the basics, and 
derive or quote some key results that are used in the main text. 


II. VARIATIONAL PRINCIPLE 


We start with a set of partial observables that constitute a {D + V)-dimensional Euclidean 
configuration space C. A point q in the configuration space represents simultaneous measurement 
of all partial observables, e.g., q = {x, 4>). To establish a physical theory, one has to specify corre¬ 
spondence between the partial observables, and physical measuring devices, such as clock, rulers, 
or instruments measuring components of the field. In this article we take such correspondence 
for granted, as we are only concerned with the mathematical aspects of the theory. 

Denote by D the dimensionality of motions, i.e., submanifolds 7 of the configurations space C. 
With D = 1 we can study particle mechanics, with D = 2 we can do string theory, or field theory 
in two dimensions, and so on. We shall not deal with systems with gauge invariance, for which 
the mathematical motion (the surface in C) has higher dimensionality than the actual physical 
motion (the trajectory). 

Tangent space of 7 at point q is spanned by D linearly independent vectors oi,..., od, which 
are conveniently combined into a grade-ZI multivector oi A ... A qd- Normalized version of the 
latter is called the unit pseudoscalar of and it is denoted by I-y. In the terminology used in 
Ref. [ini Ch. 6 ], the function I^{q) represents a Zl-dimensional distribution on C, with 7 being 
its integral manifold. 

Fundamental for the following formulation of dynamics is the concept of generalized momen¬ 
tum, which is a grade-D multivector, denoted by P, defined at each point of 7 (see Fig. [I). 
It serves as a quantity conjugated to thus generalizing the canonical momentum of particle 
mechanics. 

The last ingredient is the Hamiltonian H{q,P), which is supposed to be scalar-valued. (Gen¬ 
eralization to the case of multicomponent H is straightforward.) 

The variational principle can now be stated as follows (cf. [U Ch. 3.3.2]): 
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Variational principle. A surface 7 ci with boundary d^c\ is a physical motion, if the couple 
( 7 cb^ci) extremizes the (action) functional 

A[^,P]= [ P{q)-driq) (1) 

in the class of pairs ( 7 , P), for which dj = djci, and P defined along 7 satisfies 

H{q,P{q))=0 yqGj. (2) 


(The subscript “cl” stands for “classical” motion, or trajectory, which we sometimes use instead 
of the expression “physical” motion.) 

The integral in 0 is defined in ( |A43[ ) without having recourse to a parametrization of the 
surface 7 . Of course, if desired, the oriented surface element dr can be cast, using arbitrary 
coordinates on 7 , as 


The int egran d in 0 may be rephrased as a differential form 9 = pj,^,,,j^dq^^ A ... A dq^^ (ggg 
Formula ( |A^ in Appendix [A|). However, the merit of geometric calculus consists, in fact, in 
splitting of the differential form into two parts, dT and P, where P is able to enter into functions 
such as H{q, P). 

Finally, let us note that in (TJ Ch. 3.3.2] the action is an integral of 6 over the submanifolds 
of the bundle of D-forms over C. Since we hesitate to work in spaces mixing points q and 
multivectors P, we rather operate with surfaces in C, on which the momentum field is defined. 


III. CANONICAL EQUATIONS OF MOTION 

We will now derive the equations of motion that follow from the variational principle of Sec¬ 
tion |n[ For this purpose, we incorporate the Hamiltonian constraint ([^ into the action 0 by 
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means of a scalar Lagrange multiplier A. The augmented action is the functional 

Ah,P,X]= f [P{q)-dr{q)-Xiq)Hiq,P{q))], 

J y 


(4) 


where A is in fact an infinitesimal quantity with magnitude comparable to |(ir| - the magnitude 
of dr (see definition ( [All )). 

Varied action A['y',P', A'] is the integral taken over a new surface 7 ', and featuring new func¬ 
tions P' and A', which are defined along 7 ' (see Fig.[^. Let /, where 


/(?) = q + Sqiq), 


(5) 


be the infinitesimal diffeomorphism mapping between 7 and 7 ', i.e., 7 ' = {q' = f{q) \ q G 7 }, and 
denote by 


SP{q) = P'ifiq)) - P{q) and 5X{q) = X'{f{q)) - X{q) 

the variations of momentum and Lagrange multiplier, respectively. 

Variation of the action 5A = A[i, P', A'] - A[ 7 , P, A], is given by 


( 6 ) 


<5Al= / [P'{f{q))-f_{dT{q)-q)-X'{f{q))H{f{q),P'{f{q)))]- [P{q) ■ dnq)-X[q)H{q, P{q))] , 

J y J -y 

(7) 


where we have employed the integral substitution theorem (A44) to transform the integral over 
7 ' into an integral over 7 . For the infinitesimal diffeomorphism /, the outermorphism mapping 
/ that specifies the transformation rule for multivectors, is given by Formula (A36). Therefore, 
up to first order in Sq, SP and SX, we find 

6A= j [(P -f 6P) ■ {dr + (dr • dq) A Sq) -{X + SX)H{q + Sq,P + 6P)-P- dF + XH{q, P)] 

J y 


—SXH{q, P) + SP ■ (dF — XdpH{q, P)) — XSq ■ dqH{q, P) + P ■ ((dF • dq) A Sq) 


( 8 ) 


where the vector derivative dq, and the multivector derivative dp are defined in (^A23), and (A38), 


respectively. The “overdot” notation is used here to indicate the scope of the differential operator 
dq, and has nothing to do with time derivative. Without an overdot, any differential operator is 
supposed to act on functions that stand to its right. 

The last term in (|^ can be recast with a help of the Fundamental theorem of geometric calculus 

(ira. 


f P ■ ((dF ■dq)XSq)= j P • (dE ^Sq)- f P • ((dF • dq) A Sq), 

J y J Q-y J 'y 


(9) 


where dE is the oriented volume element of the boundary d^. Now, the first term on the right- 
hand side vanishes, since we assume that 7 and 7 ' have common boundary. 

For D — 1, dr ■ dq \s algebraically a scalar, and so the integrand in the second term is readily 
reshuffled. 


P • ((dF • dq) A Sq) = Sq ■ (dF • dqP). 


( 10 ) 


(Mind the priority of the inner product and the outer product “A” before the geometric 
product, which is denoted by an empty symbol.) 
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For D > 1, we can use identities ( |A10d ) and ( |A10a ) to find 

p ■ ((dr • dg) A Sq) = {p ■ (dr • dg)) ■ Sq = (-i)^-Mg • ((dr • dg) ■ p). 


( 11 ) 


The two cases have to be treated separately due to the definition ( |A 6 [ ) of the inner product. 
After these rearrangements we arrive at our final expression for the variation of the action, 


5A'. 


-5X H{q, P)+6P- (dr - A dpH{q, P)) +Sq-{ (-l)^(dr ■dg)-P-\ dgH{q, P) 


( 12 ) 

which holds for D > 1, while the case D = I is obtained simply by replacing (dr • dg) ■ P with 
dr • dgP. The requirement that the variation of the action be zero for all SP, 6q and 6X yields 
the following 

Canonical equations of motion. Physical motions yd are obtained by solving the system of 
equations 


XdpH{q,P) = dr, 


{-l)^XdgHiq,P) = J'^^ « 


dr • dgP for Z? = 1 

(dr ■dg)-P for D>1, 


H{q,P)=0. 


(13a) 

(13b) 

(13c) 


(We use the adjective “canonical”, because these equations generalize, as we shall see in Ex¬ 
ample V A Hamilton’s canonical equations of motion of classical mechanics.) 


The first canonical equation (13a) furnishes a relation between the momentum P, and the 


tangent planes of 7 represented by the oriented surface element dT. It asserts that the multivector 
derivative dpH, which is a grade-Z? multivector, is proportional to dr, with the proportionality 
constant equal to A. Note that one can always normalize dr and A by the magnitude |dr| to free 
Eqs. (13) from infinitesimal quantities. 


The second canonical equation (13b) describes how the momentum multivector P changes as 


it slides along the surface 7 . It is important to note that P is differentiated effectively only in 
the directions parallel to 7 , as a consequence of the inner product between the surface element 
dr, and the vector derivative dg. Moreover, the “overdot” on the left-hand side assures that only 
explicit dependence of ZZ on g i s be ing differentiated, not the dependence through P{q). 

The last canonical equation (13c) is the Hamiltonian constraint (§. Let us remark that had 


we started with several constraints Hj(q,P) = 0 in the variational principle, we would have 
introduced the corresponding number of Lagrange multipliers Xj, and the canonical equations 
would contain the terms instead of AZZ. 

In Appendix we provide a component form of the canonical equations to make them more 
accessible for a reader who is not sufficiently familiar with the formalism of geometric algebra. 
Note, however, that this step is not necessary in order to make practical calculations, as will be 
illustrated in the examples of Sec. |V] 


IV. LOCAL HAMILTON-JACOBI THEORY 


One possible method to approach the canonical equations (13) is the following. Suppose P{q) 
is given, that obeys the Hamiltonian constraint 


H{q,P{q)) = 0 


( 14 ) 
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on some open subset of the configuration space C. By differentiation, we obtain, according to 
the chain rule (A41), 


dgH{q,P{q)) + dqP{q) ■ dpH{q,P{q)) = 0, 


and using the first canonical equation (13a), we find 


XdqH{q,P{q)) = -dqP{q) -dV. 


(15) 


(16) 


But the right-hand side may be rearranged by means of identity (A8) for H = 1, or (AlOal and 
(AlOg) for H > 1, with the result 


XdqH{q,P{q)) = 


dT ■ {dq A P{q)) — dT ■ dqP{q) for Zl = 1 

(-l)^-idr- {dqXP{q))+{-l)^{dP-dq)-P{q) for £» > 1. 

Therefore, we observe that if 


(17) 


dq A P{q) = 0, 


(18) 


the second canonical equation (13b) is automatically fulfilled. Momentum satisfying this condi¬ 
tion can be expressed, at leas t locally, as P{q) = dq A S{q), where 5” is a multivector of grade 


D — 1 (consider relation (A27)). Canonical equations (13) are then reduced to two equations, 

(19) 


XdpH{q,dqXS) = dT 
and the local Hamilton-Jacobi equation 

H{q, dqAS)=0 


( 20 ) 


(see Appendix for a component form of this equation). If we succeed in finding a solution of 
Eq. (pOl), we can plug it into Eq. (19), which then defines a distribution of tangent planes of 


a classical motion surface. This distribution is integrable only if certain conditions are satisfied 
(see [ini Ch. 6.1]). 

In addition, if we find a whole family of solution S{q',a), parametrized by a continuous pa¬ 
rameter a, then differentiating Eq. (20) with respect to a, and substituting Eq. yields 


0 = Xd^H{q,dq AS) = Xda{dq A 5) • dpH{q,dq A S) = dT ■ {dq A (daS)). 
Now, for D = 1, the Hamilton-Jacobi function S is scalar-valued, and we have 
dr • dq{daS) =0 ^ daS{q; a) = P Vg G yd, 


( 21 ) 


( 22 ) 


for some constant /3, meaning that the quantity daS{q; a) is conserved along physical motion. If 
one finds N such parameters (recall that the dimension of the configuration space is A^ -I- 1), the 
physical motion yd can be given implicitly by a set of constraints between partial observables, 

da^S{q-ai,...,aN) = di 


) Q^iv) = Pn- 


(23) 


Of course, we assume that the N constraints are independent, i.e., that the gradients 
dq{da^S ),..., dq{da^S) are at each point linearly independent. In Example VC we will illustrate 
the Hamilton-Jacobi method with the case of a relativistic particle. 
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When D > 1, we can use identity (AlOd) and theorem (A45) to cast Eq. (21) as 


(dr • dq) ■ (daS) = 0 ^ f {dr-dq)- (a„S') = / dS • {d^S) = 0, 

dlc\ 


(24) 


where 7 ci is an arbitrary H-dimensional subset of 7 ci (a “patch” on 7 ci). Therefore, in the 
multidimensional case, the conservation law (22) is replaced with a certain continuity equation. 


One remark is in order before closing this section. In classical particle mechanics, one of the 
solutions of the Hamilton-Jacobi equation is the action along classical trajectory, regarded as 
a function of one of the endpoints. In field theory, the classical action may be viewed as a 
functional of the boundary d"fc\- Some authors (e.g. [U Ch. 3.3.4]) have therefore considered a 
variational differential equation that describes how the classical action changes under variations 
of the boundary, naming it also “Hamilton-Jacobi equation”. Note that Eq. (20) is substantially 


different from this kind of approach, for it contains only partial, not variational, derivatives. 
That is why we call it “local Hamilton-Jacobi equation”. Local Hamilton-Jacobi theory is also 
treated, e.g., in Refs. [5] and [3]. 


V. EXAMPLES 

In the following examples, we illustrate the general theory by specifying a concrete form of the 
Hamiltonian H{q, P). Reader’s familiarity with the techniques of geometric algebra and calculus 
on the level of Appendix fA] is assumed. 


A. Non-relativistic Hamiltonian mechanics 

Let us consider U = I, and choose a constant unit vector e* in the configuration space 
C ~ Arbitrary point q can be decomposed as q = t + x, where t is parallel to et, 

while X is perpendicular to et (see Fig.[^. 

Define the Hamiltonian as follows: 


H{q,P)=P-et + Ho{q,P), (25) 

where Ct ■ dpHo = 0. We shall identify Hq as the non-relativistic Hamiltonian of a mechanical 
system. 














dotting” the first canonical equation (13a I with e*, we find 


et-dV = \efdpH = (26) 

This means that the tangent vector of a physical motion yd has nonvanishing component parallel 
to et, and we can therefore present 7 ci as 

7ci = {<7 = g{i) = t + x(t) 1 1 G spanjej ~ K}. (27) 


The line element on the t-axis is related to the line element on yd via the differential mapping 

(IMsl), 


Eq. (26) then implies 


dr = g{dt) = dt + dt ■ dtx. 


(28) 


X = et ■ dT = et ■ dt = etdt ^ dt = Xet ^ dT = X g{et) = X{et + e* • dtx), (29) 


which allows us to eliminate A from the equations of motion. 

Let us denote p{t) = P{g{t)), and observe that the first canonical equation (13a) assumes the 
form 


et ■ dtX = dpHo{q,p). 


(30) 


As concerns the second canonical equation ( 13b[ ), we realize that by the chain rule for differ¬ 
entiation 


et • dtp{t) = et ■ dtP{g{t)) = g{et) ■ dqP{q)\g=g(^t) = jdT ■ dqP{q)\q^g^t), 


and so we arrive at 


et • dtp = -dqH{q,p) = -dqHo{q,p). 


Projecting Equation (32) onto ex, an arbitrary vector perpendicular to et, we find 

et * dt ex ' p ex * dqPQ {q, p), 


while projecting onto et yields 


et- dtefP = -et ■ dqHo{q,p). 


Using finally the Hamiltonian constraint, Eq. (13c), with H given by (25), Eq. (34) implies 

et • dtHo{q{t),p{t)) = et ■ 9gido(<7,p(i))lg=g(f)- 


(31) 

(32) 

(33) 

(34) 

(35) 


It is now easy to realize that Eqs. (30) and (33) represent Hamilton’s canonical equations 
of motion for a non-relativistic system with non-relativistic Hamiltonian Hq, while Eq. (35) 
expresses the law of conservation of the total energy Hq, if et ■ dqHo = 0, i.e., if Ho{q,p) does not 
depend explicitly on time. Intuitively, the Hamilton’s canonical equations follow from variations 
of the trajectory 7 in the x-space, while the energy conservation is a result of variations in the 
t-space (see Fig. [^. 

Hamilton-Jacobi equation (20) for a scalar function S{q) reads 

H{q,dqS) = et ■ dqS + Ho{q,dqS) = 0, (36) 


and hence reproduces the standard Hamilton-Jacobi equation of classical mechanics. 
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B. Scalar field theory 


In this example we will show that the formalism based on Hamiltonian constraint can accom¬ 
modate the theory of real one-component scalar field > 5 ^( 3 :), defined on a U-dimensional Euclidean 
spacetime by the Lagrangian 

C{ct>,d,ct>) = ]^{dAf-V{ct>). (37) 

For this purpose, let us assume D > 1, choose a unit Z3-blade Ix in a (I? -I- l)-dimensional 
configuration space C, and define the Hamiltonian 


H{q,P)=P-Ix+HDw{q,P), 


(38) 


where Ix ■ dpHow = 0. The blade Ix defines a splitting of the configuration space C into a 
£)-dimensional spacetime, spanned by an orthonormal set of vectors {ei,... , 60 }, Ix = ei ■ ■ - Cd, 
and the field space, which is its one-dimensional orthogonal complement, represented by a unit 
vector By (see Fig. [^. 

In analogy with Example V A let us take the inner product of the first canonical equation 
(13aI with Ix, the reversion of Ix- 


Ix ■ dr = XIx ■ dpH = XIx • /a; = A yf 0. 

This implies for any vector a tangent to 7 ci (i.e., a A dT = 0) that 

0 ^ Ix ■ dr = Ix ■ [a A {a~^ ■ dT)] = {Ix ■ a) ■ (a~^ ■ dT) ^ Ix ■ 0. 


(39) 


(40) 


This means that no tangent vector is perpendicular to Ix, and we can therefore assume that a 
physical motion yd is presented as 

7ci = {q = g{x) =x + y{x) I X e span{ei,..., ep,} ~ K^}. 


(41) 
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The surface element dT is related to the infinitesimal element of the cc-space dX = \dX\Ix via 
the outermorphism mapping induced by the function g, dT = g{dX). It is not hard to show, in 
analogy with Equations (A35l and (A36), that for any r-vector from the a;-space, 


g{Aj.) — Aj. + (Aj. • dx) A y, 


where dx = ' dq. Eq. (391 then supplies the relation 

\ = Ix ■ dV = IxdX => dX = XIx => dT = Xg{Ix) = X[lx + {Ix ■ dx) A y). 


The first canonical equation (13a) reduces to 

{Ix • dx) Ay = dpHow- 


(42) 


(43) 


(44) 


As concerns the second canonical equation, we observe that by the chain rule for differentiation, 
Cj • dxP{g{x)) = g{e.j) ■ dqP{q)\q—g(^x) = ■ 9{dq)P{<])\q=g(x) ^ dx = g{dq), (45) 


where g is the adjoint of the linear mapping g (see Eq. (A30l). With a usual abuse of nota tion, 
we will regard P as a function of x, P{x) = P{g{x)), from now on. Using the rules ( A33[ ), we 
can cast the left-hand-side of the second canonical equation (13b) as g{Ix ■ dx) ■ P, which, by 
Formula (42), yields 

{Ix - dx) ■ P + {[(4 • dx) ■dx]Ay}-P= {-l)^dqH = {-l)^dqHDw, (46) 

where the accent has the same role as the dot in that it indicates the scope of differentiation of 
the corresponding diffe rent ial operator. 

Now, “dotting” Eq. (44) with Ej, where Ej = Ixejey, so that Ej -4 = 0, we find, after little 
(geometric) algebra 


Cj ■ dx Cy ■ y = Ej ■ dpHow: 


(47) 


while Eq. (46) multiplied by e„ reads 


Cj * dxEj 


P = —e-y ■ dqHpw 


(48) 


(Summation from 1 to Z? over a repeated indices is implied here and below.) Notice that Eqs. 
(47) and (48) are identical to the De Bonder-Weyl equations of motion, while Hpw is identified 
with the De Bonder-Weyl Hamiltonian [Mi- 
Taking specifically 


1 ^ 


i=i 


where cj) = Cy ■ y, we obtain from Eq. (47) the relation 


D 


Cj ■ dx(j) = '^{Ej ■ Ek){P ■ Ek) = P ■ Ej, 


k=l 


(49) 


(50) 


which substituted in Eq. (48) yields 


dlct> = -d^V{ct>). 


( 51 ) 
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This is the equation of motion of a scalar field described by the Lagrangian (37). 

In addition, taking inner product of Eq. (46 1 and a spacetime vector provides an interesting 
relation analogous to the energy conservation, Eq. (35), of non-relativistic mechanics. We have 

[cj A (4 • d^)] ■ P + {Cj A [(4 • dj:) P = • dqHow, 


and introducing the resolution P = I^Ix ■ P + EkEk ■ P we get 

[cj A (4 • dx)] ■Ixlx-P + {ej A [(4 • dx) ■ dx] A y} ■ E^ E^ ■ P = (-1)^6^ • dqHow- 
After some algebra we obtain (recall that Ej = eyejlx) 

- Cj ■ dxP ■ 4 + {-l)^{ej A [(4 • dx) ■ {dxft))]} ■ {ekTx)Ek ■ P = Cj ■ dqHow- 
The second term on the left-hand side is simplified as follows: 

(-l)^{ej A [(4 • dx) ■ {dxc!))]} ■ {ekIx)Ek ■ P = 

= {ej A Cfe A [(4 • dx) ■ {.dx(t>)\} -IxEk- P 
= {cj A Cfc) ■ {[Ix- {dx A {dx(j)))] ■ Ix}Ek ■ P 
= {ej A Cfc) • {dx A {dx(l)))Ek ■ P 
= {cj A Cfc) • {dx A {dx^))Ek ■ P 

= -ej ■ dx[{ek ■ dx(j)) {Ek ■ P)] + ek ■ dx[{ej ■ dxcj)) {Ek ■ P)], 

where the second last equality is a consequence of commutativity of partial derivatives. 
Using in Eq. (54) the Hamiltonian constraint (13c), with H given by (38), we arrive at 

Cfc ■ dx[^jkEDW ^jk{ei * dx({'){Ei ■ P) -\~ {ej ■ dx4^){Ek * T*)] — ej * dqE[j^\^. 


(52) 

(53) 

(54) 


(55) 


(56) 


For H£)\Y defined in (49) we can eliminate the momenta, owing to relation (50), and obtain 

ek ' dx\ ^jkE{4^: dx4^) {ej • dx4^){ek ■ dx4^)\ — ej * dqH — 0 . (^^) 

The expression in square brackets is nothing but the canonical energy-momentum tensor Pjk of 
the scalar field with Lagrangian (37), so we finally arrive at the continuity equation 

ek ■ dxTjk = 0. (58) 

Note that the equation of motion for the scalar field, Eq. (51), and the continuity equation 
for its energy-momentum tensor have common origin in the canonical equations of motion (13). 
Intuitively, the equation of motion is a result of variations of the surface 7 in the y-direction, 
while the continuity equation follows from variations in the x-plane (see Fig. |^. 

Finally, let us say a few words about the Hamilton-Jacobi theory. Equation ( [20| with Hamil¬ 
tonian (38) reads 


4 • {dq A S) + Plow {q,dq A S) =0, 


(59) 


where S{q) is a multivector of grade D — 1. For concreteness, consider Hq given by Eq. (49), and 
assume S is parallel to 4- Defining the vector s{q) = S{q) ■ 4j also parallel to 4) we observe 
that 


4 • {dq A S) = dq - s, 

Ej ■ {dq A S) = ey ■ dq ej ■ s. 


(60) 


Eq. (59) then takes the form (note that dj,s = ey ■ dqs) 

dq-s + ]^ (90s)^ + V{4>) = 0, (61) 

which coincides with the field-theoretic Hamilton-Jacobi equation derived formerly by Weyl 
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C. String theory 


Probably the simplest nontrivial Hamiltonian to consider is 

H = \{\P\^ - 


where A > 0 is a scalar constant, and |P| is the magnitude of P (see definition (All)). 
According to Formula (A40), the first canonical equation (13a I takes the form 

dr = AP, 


( 62 ) 


(63) 


which substituted into the Hamiltonian constraint (13c I fixes the absolute value of the Lagrange 
multiplier A, 


|dr| = |A|A. 


(64) 


by A, and using Eq. (64) we find 


Furthermore, substituting Eq. (63) into the second canonical equation of motion (13b), dividing 


■ dq = 0 {D = 1), 
{I^ ■dq)-I^ = 0 {D> 1), 


(65) 


where I-y = dr/|dr| is the unit pseudoscalar of the surface 7 . This equation has a simple 
geometric interpretation. It entails vanishing of the mean curvature of 7 , or of its generalization, 
the spur vector (see Ref. O Ch. 4-4]). 

Eqs. (63) and (64) allow us to rewrite the action 0 in terms of dF, 


f P dT= f i|dr|2 = ±A f |dr|, 

J ■y J -y ^ J 'y 


( 66 ) 


where “±” is the sign of A, and |dr| = VdF • dF. This is the Euclidean Nambu-Goto action of 
the string theory m- It is proportional to the volume of the worldsheet 7 , with A playing the 
role of string tension (the speed of light is set to unity). Extremals of this action, i.e., solutions of 
Eq. ( |65| ), minimize the volume, and so are called minimal surfaces in mathematical literature [ 8 ]. 
It is worthwhile to mention that the Nambu-Goto string can be formulated also within the De 
Bonder-Weyl Hamiltonian theory (see Ref. [T^], which uses the language of differential forms). 

In this example the Hamilton-Jacobi equation (20) takes a particularly compact form (cf. 
Refs. [13] and [^ Ch. 7]) 


\dq AS\=A. 


(67) 


Erom now on, let us focus on the case D = 1, which describes the relativistic particle in the 
Euclidean spacetime. We will present two methods for finding the physical motions. 

First, suppose that two points, qo and q, lie on 7 ci, multipl y Eq . (65) by |dr|, and integrate 
along 7 ci from qo to q. The Fundamental theorem of calculus (A45) implies that 

I^{q)-I^{qo)=0, ( 68 ) 

i.e., Ij is constant along a physical motion, and are therefore straight lines in C, 


lc\ = {q = VT P qo\T & 


(69) 


where qo € C and v is arbitrary constant vector. 
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Second method utilizes a family of solutions of the Hamilton-Jacobi equation (67), for example, 

Siq;qo) = 9o\- (70) 


According to Formula (22), derivative of S with respect to the parameters qo yields conserved 
quantities 


5,0 5 = 


-A 


g- go 
|g- gol’ 


Physical motion are then obtained readily, 


7ci = 


g-go 

|g-go| 



where v is an arbitrary constant unit vector. 


(71) 


(72) 


VI. CONCLUSION AND OUTLOOK 


In this article we elaborated on the formulation of classical field theory presented in [1] Ch. 3], 
which is based on the notion of partial observables, and on the Hamiltonian constraint. The latter 
is a function of configuration-space point and the generalized multivector-valued momentum. 
Starting from the variational principle of Section |ll] we derived canonical equations of motion 
( [T^ . We also deduced local Hamilton-Jacobi equation (20), which can be a useful tool to find 
the physical motions. 

With three ensuing examples we showed how non-relativistic mechanics, scalar field theory, 
and string theory can be described in one unifying framework by appropriately selecting the 
Hamiltonian constraint. In particular, we noticed that equations of motion, and the continuity 
equation for the energy-momentum tensor (which reduces to the energy conservation equation in 
the case of non-relativistic mechanics) are in fact of the same origin in the Hamiltonian constraint 
formalism. Therefore, this formalism may be of interest even for theories that do not assume 
symmetry between time and space, or spacetime and fields. 

Although we restricted our attention to Euclidean space, an extension of the formalism to 
pseudo-Euclidean spaces should be relatively straightforward [H Ch. 1-5]. Then, one has to 
mind, and keep track of, possible sign differences between the reversion and the inversion of unit 
multivectors. Eor example, in general, ^ I^. 

Hamiltonian formalism is especially important when it comes to quantization. In particle 
mechanics, momentum is promoted to a differential operator, and the Schrodinger equation is 
postulated. What is the quantum operator corresponding to the multivector-valued generalized 
momentum of the Hamiltonian constraint approach? And what does the Schrodinger equation 
look like, once we know the classical Hamilton-Jacobi equation (20). Although these questions 
have not been addressed in general, let us note that there have been studies of quantization 
in the De Bonder-Weyl Hamiltonian theory, where the quantization of momenta is based on 
generalized Poisson brackets, and a field-theoretic generalization of the Schrodinger equation 
is proposed that features a Clifford-valued wave function, and reduces to the De Bonder-Weyl 
Hamilton-Jacobi equation in the classical limit mm- 
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Appendix A: Geometric algebra and calculus 

We give a brief introduction into the formalism of geometric algebra and calculus in a way 
that respect the requirements of this article. For a thorough and rigorous treatment, the reader 
is advised to consult monograph [^, which we shall frequently quote. Gomplementary to this is 
the textbook |15j . which provides, apart from a pedagogical introduction into the mathematical 
formalism, many diverse physical applications. 


1. Geometric algebra 


Let us start with an n-dimensional real vector space V, and define the geometric product of 
vectors by the following axioms: 

Va, b,c : 1) a(bc) = (ab)c = abc 

2 ) a(b -I- c) = ab + ac 

3) > 0 for nonzero a. (Al) 


This product induces an associative algebra over vector space V — the geometric (or Clifford) 
algebra G{V). Frequently, the term “Clifford algebra” can be encountered in literature. Never¬ 
theless, we prefer the name “geometric algebra” used originally by Clifford, and advocated by 
Hestenes [9] to emphasize its geometric interpretation. Let us note that although it is possible 
to represent vectors in the algebra by Dirac gamma matrices, it is in fact not very useful, since 
it provides no insight int o the properties of the algebra, neither it simplifies any calculations. 

The last axiom in (Al) has far-reaching consequences. Expanding (a -I- b)^ = + b^ + ab + ba, 

we observe that 


a • b := 2 


(A2) 


is a scalar. The remaining part of the geometric product, 

a Ab := — ba), 


(A3) 


is a bivector. The geometric product of two vectors is therefore decomposed into two parts: 
symmetric non-associative inner product (A21, and antisymmetric associative outer product (A31. 
The scalar a • 6 is identified with the scalar product of vectors a and b. Positive definiteness of 


this scalar product follows from the strict inequality in the last axiom in (All. One could also 


consider indefinite quadratic forms, but we do not deal with them in this text. 

The bivector a Ab represents an oriented parallelogram spanned by the two vectors. In fact, 
it represents a whole equivalence class of parallelograms, since, e.g., aAb = aA{b-\- Xa) for 
arbitrary scalar A, as follows from the antisymmetry of the “A”-product. 

Successive multiplication of vectors generates the entire geometric algebra. General elements 
are called multivectors. They decompose into a sum of terms with different grade. An element 
of the algebra has grade r (it is an r-vector), if it can be written as an exterior product of r 
vectors, in which case it is called r-blade, or as a linear combination of such terms. Because 
Ar = oi A ... A Or vanishes if and only if the vectors Oj are linearly dependent, it represents 
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an object with dimensionality r, namely, the parallelotope spanned by the vectors Oj, whose 
orientation is specified by the order in which the vectors appear in the exterior product. By the 
Gram-Schmidt process [HI Ch. 1-3] it can be shown that every blade = Oi A ... A a^. can 
be written as a scalar multiple of a geometric product of orthonormal vectors: Ar = a ei... Cr, 
where ej ■ Ck = Sjk- 

General multivector is a linear combination of terms with increasing grade: 

1 , Cj , (j ^ 5 ■ ■ ■ ; €i . . . 

where j,k,... = 1,... ,n. Geometric algebra G{V) is therefore the linear span of these terms. A 
useful convention states that scalars are denoted by Greek letters a, /3,..., vectors are denoted by 
lower case Latin letters a,b,.. and other multivectors are denoted by capital letters A,B,.... 
Note, however, that in physical applications there can appear exceptions to these rules due to 
conventions used in physics. For example, the Hamiltonian in ([^ is denoted by H although it is 
a scalar function. 

When a multivector has some definite grade r, this fact is indicated by a subscript, e.g., we 
write A,.. The geometric product between two such multivectors decomposes into a sum of terms 
with specific grade, 


A^Bg — |^_g| -t- |,._g|_|_2 “b ... ^Aj-Bg)j-j^g^ 


(A5) 


where the symbols {M)j. denotes the projection onto the grade-r component of M. Projection 
onto the scalar part is abbreviated (M) = {AI)q. Lowest- and highest-grade terms in the series 
(A5) are designated by 


Ar ■ Bg 


if r, s > 0 

0 if r = 0 or s = 0, 


Ar A Bg — f^ArB g^ r+s-! 


(A6) 


which are again called interior and exterior product, respectively. The inner product is defined to 
be zero if either of the multivectors is a scalar. If one of the members is a vector, then relations 
(|A6|) reduce to 


a ■ Ar = -(aAr — (—lyArU), 
a A Ar = “ {^aAr -t- (— 1) . 


(A7) 


Moreover, expressions such as a • A can be expanded in primitive terms using the identity (see 
[i Ch. 1-1]) 


a ■ (oi A ... A ar) = ^^(“1)'^ • Oj Oi A ... A d^- A ... A a^, (AS) 

where the check marks vectors that are dropped out from the expression. 

To avoid overload of brackets, we have adopted the standard convention that inner and outer 
products have priority before the geometric product. For example, 

a • A 6 A c = (a • A)(6 A c). (A9) 

To every blade A,. = ai A ... A corresponds a subspace of V spanned by the vectors aj. A 
generic vector a € V belongs to this subspace if and only if it is a linear combination of aj’s, i.e.. 
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if and only if a A = 0, as follows from antisymmetry of the outer product. Alternatively, the 
condition a A = 0 can be restated as a - Ar = aA^. Vectors in the orthogonal complement are 
characterized by the requirement a ■ A^. = 0 inferred from expansion (A8). 

We quote several useful identities derived in [21 Ch. 1-1] that enable efficient manipulations: 


Ar ■ Bs = (—I)’'*-'* ■ Ar for r < s, (AlOa) 

ArABs = {-lY^B,AAr (AlOb) 

Ar ■ {Bs ■ Ct) = {Ar A Bs) ■ Ct ioT r + s <t and r, s > 0, (AlOc) 

{Ct ■ Bs) ■ Ar = Cf {Bs A Ar) for r + s <t and r, s > 0, (AlOd) 

Ar ■ {Bs ■ Ct) = {Ar ■ Bs) ■ Ct foT r + t < s, (AlOe) 

a ■ {Ar A Bs) = {a ■ Ar) A Bs + {—lyAr A (a • Bs), (AlOf) 

a A {Ar ■ Bs) = {a ■ Ar) ■ Bs + {—lyAr ■ {a A Bs) for s > r > 1. (AlOg) 


To every multivector A is associated a scalar magnitude |A| by 


|Ap = (AA), 


(All) 


where ” . ” is the operation of reversion, defined trivially on vectors, a = a, and extended to 
G{V) by linearity and the requirement 


AB = BA. (A12) 

For a blade Ar = ai A ... A Or, |A| is indeed the volume of the parallelotope spanned by the 
vectors aj. Also note that 


Ar = (-1)A’'-i)/2a^, 


and since ArAr 


{ArAr), the blade Ar has an inverse. 


A 


-1 

r 




(A13) 


(AM) 


For unit blades, the inverse is equal to the reverse: A,, ^ = Ar. 

Every grade-r multivector gives rise to a scalar-valued function a of r vector variables, 

a{bi,... ,br) = Ar ■ {bi A ... A br), (A15) 

which is linear in each argument, and changes sign whenever two vectors are exchanged, i.e., a 
is an alternating form. In fact, every alternating form can be represented by some multivector 
Ar in this way, and operations on forms can be naturally expressed in terms of operations on 
the corresponding multivectors 0 Ch. 1-4]. 

If Ar is a blade, then the inner product of two r-blades, Ar = Or A.. .Aai and Br = bi A.. .Abr, 
can be expressed as the determinant of the matrix of scalar products Oj ■ bk. 


Ar ■ Br = det(aj • bk). 


In addition, using the expansion (A8) we derive 


n 

{Ar ■ a) ■ {b ■ Br) = ^ {-iy^"'{a ■ ai){b ■ bm)^i^OT{aj ■ bk\l,m), 

l,m—l 


(A16) 


(A17) 
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where 


minor(aj • m) = (a^ A ... A a; A ... A oi) • (5i A ... A A ... A &r) (^18) 

denotes the {l,m) minor of the matrix aj ■ bk, i.e., the determinant of the latter matrix with the 
Zth row and mth column erased. 

Let us choose an orthonormal basis of the vector space V. Any vector a can be 

expanded as 

n 

a = '^{a-ej)ej, (A19) 

1=1 

where a ■ Bj are the components of a with respect this basis. Moreover, any grade-r multivector 
Ar can be expressed as a sum 

Ar = ^ {Ar ■ ej)ej = ^ {Ar ■ ej)ej (A20) 

\J\=r \J\=r 


over a ll ordered sets of indices J = (ji, ... ,jr), ji < ■ • ■ < jr, where ej = ... 

(A20) can be proved with a help of the formula (see Eq. (3.14) in [9l Ch. 1-3]) 


Relation 


ck -ej = {ckr ■ ■ ■ efej • (cji ... e^-J = = Sf. 


(A21) 


2. Geometric calculus: Differentiation 

We will now move towards the geometric calculus, that is, the theory of differentiation and 
integration developed by D. Hestenes [5], which takes advantage of the rich algebraic structure 
of geometric algebra. 

Our setting involves a Euclidean vector space V, corresponding to the configuration space C 
of partial observables, which holds at every point q € V a, copy of geometric algebra G{V). Take 
function F{q) with values in G{V), and a vector a. The derivative of F in the direction a is 
defined in the usual manner, 


dqF{q) ■- lim 


F{q + ea) - F{q) 


s —^0 


(A22) 


The vector derivative of F is defined with the help of an orthonormal basis of 


dqF{q) ■= ej{ej ■ dq)F{q) 


(A23) 


(summation over j is implied). The operator dq = ej{ej ■ dq) has algebraic properties of a vector, 
and hence we may separate the vector derivative into two parts. 


^qF = ^q■F + ^q^ F, (A24) 

called divergence and curl, respectively. When F is vector-valued, we recover the familiar differ¬ 
ential operators of vector calculus. For scalar F, dqF is simply the gradient. 

The vector derivative is obviously linear. To express the product rule we employ the “overdot” 
notation specifying which function in the product is being differentiated, 

dg{FG) = dgFG + d,FG. 


(A25) 
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The reason is that the vector dq need not commute with the other multivectors. Of course, 
overdots can always be eliminated by introducing a basis, for example, dqFG = ejF{ej ■ dq)G. 
Some explicit formulas for vector derivatives of elementary functions can be found in [Sj Ch. 2-1]. 

Differential forms are skew-symmetric linear functions of differential arguments dqi ,..., dqr 
that may vary from point to point. Just like in Eq. (A15), they can be expressed as 


a{dqi A ... A dq^.; q) = Ar{q) ■ {dq^ A ... A dq^.), 


(A26) 


where Ar is an r-vector function (see [9l Ch. 6-4]). Exterior derivative of a is tantamount to 
taking the curl of A^, 


da{dqi A ... A dq^_^_i-,q) = [Ar A dq) ■ {dq^ A ... A dg^_|_i). (A27) 

Let us now consider a diffeomorphism / that maps points qGVtoV, and form the directional 
derivative 


l{a-,q) = a-dqf{q). (A28) 

This gives rise to a g-dependent linear function, the differential of /, mapping vector a to a new 
vector /(a). (In standard differential geometry / is called the push-forward derived from the 
diffeomorphism /.) It is natural to extend the domain of / to general multivectors, demanding 
linearity and the outermorphism property [i Ch. 3-1] 

1{AaB)=1{A)AI{B). (A29) 

Let / denote the adjoint of the linear transform /. It fulfils, for any two vectors a and &, 

b-l{a) =J{b) ■ a =A 7{b-q) = dqf{q)-b, (A30) 

and corresponds to the pull-back of differential geometry. The adjoint is extended to an outer¬ 
morphism in the same way as the differential, 

f{AAB)=J{A)Af{B). (A3I) 

For scalar arguments we define 

/(a) = 7(a) = a. (A32) 

Although the inner product is not, in general, preserved by the differential and adjoint outer- 
morphisms, the following useful relations hold [HI Ch. 3-1] 

Ar-JiBs) =7[l{Ar) ■ Bs] for r < s, 

l{Ar) ■ Bs = f[Ar ■ 7(Bs)] for r>s. (A33) 

Let us investigate the transformations of multivectors under an infinitesimal diffeomorphism 
/(?) = 9 + bq{q). For a vector we have 

/(a) = a -\- a ■ dq5q. (A34) 

For a blade = a ei A ... A e^, where is an orthonormal basis of F, we derive 

/(Ar) = a(ei -\- ei ■ dq Sq) A ... A [cr -\- Cr ■ dq Sq) 

r 

« Ar -I- a ^ ei A ... A (cj • dq Sq) A ... A Cr 
i=i 

r 

= Ar a A .. . A Cr) • Cj] A {cj ■ dq Sq), 

i=i 


(A35) 
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where we have used backwards the expansion formula (A 81 . Now, the sum over j can be extended 
to run from 1 up to n, since e^’s are orthogonal to for j > r. Hence, applying the definition 
of vector derivative (A23) we arrive at 

f{Ar) ~ Ar + (Ar ■ dq) A 8 q. (A36) 


The Hamiltonian H{q, P) is a function of vector variable q and 11-vector variable P. In 
order to differentiate the composite function H{q, P{q)) with respect to g, we need a notion of 
differentiation with respect to the multivector variable P. 

Therefore, suppose F{P) is a multivector-valued function of a grade-1? multivector argument 
P, and A is an arbitrary H-vector. We define the A-derivative 


A ■ dpF(P) := lim 

E-fO 


F{P + eA)-F{P) 


(A37) 


and the multivector derivative 


dpF{P) := ^ ej{ej ■ dp)F{P), 

\. 1 \=D 


(A38) 


where {ej}\j\^p) is an orthonormal basis of the subspace of H-vectors of the geometric algebra 

Q{V). 

Multivector derivatives of some elementary functions are listed in 0 Ch. 2-2]. For example, 

A-dp\P\^ = A-P + P-A = 2A-P, (A39) 


which implies 


dp\P\^ =2P. 


(A40) 


The chain rule for differentiation gives 

a ■ dqH{q,P{q)) = (a • dq)H{q,P{q)) + (a ■ dq)iP{q) ■ (9p)H(g, P)|p=p(,), (A41) 

where the meaning of overdots should be evident. 


3. Geometric calculus: Integration 

Consider a H-dimensional submanifold 7 of V, whose tangent space at every point q € j is 
represented by a unit H-blade I-y{q)- I-y is called the unit pseudoscalar 0/7 [3 Ch. 4-1], and 
it defines also the orientation of 7. Recall that a vector a(q) is tangent to 7 if and only if 

Directed integral of a multivector-valued function F{q) over the manifold 7 is defined in an 
intuitive way, by approximating 7 with a chain of simplices Ar(g) (see [T^, [3 Ch. 7], or [13 
Ch. 6.4]), and taking the Riemann sum, 


/ dr{q)F{q) := lim VAr(gi)F(gi 

I n—>-00 ^ 


(A42) 


Here dT = jdrjJ.y is the oriented surface element of 7 . Undirected integral of F is equal to the 
directed integral of I~^F. 
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In a more general setting, we can define the directed integral 

« n 

/ L{dT{q);q) ■= l™ '^LiAT{q^);q^), 

»_i 


(A43) 


i=l 


where L(Ad] q) is a multivector-valued function, linear in the U-vector argument Ajj. Obviously, 
the integral (A43) reduces to (A42) if we choose L{AD]q) = ADF{q). 

Suppose we are given two surfaces 7 and 7 ', where 7 ' = {q' = f{q) | 9 G 7 } for a certain 
diffeomorphism /. The unit pseudoscalar is mapped by the outermorphism / onto 7 ', where 
it is proportional to We shall assume that the proportionality constant is positive, i.e., that 
the surfaces 7 and 7 ' have common orientation. 

Oriented surface elements are related by dT' = f{dT). To compare integrals along the two 
surfaces, we can change the integration variables (see also O Ch. 7-5]), 


[ L{dT'{q');q')= [ L(/(dr( 9 ); g);/(q)). 
J y' J 'y 


(A44) 


Finally, let us quote the Fundamental theorem of geometric calculus [g Ch. 7-3], [m Ch. 6.5], 
or the generalized Stokes theorem, which relates integral over a surface 7 and integral over its 
boundary dj: 



L{dT,-, q) 


L{dT ■d;q), 


(A45) 


where here the first argument of L is a multivector of grade D — 1, and dS is the surface element 
of the boundary dj. 


Appendix B: Component form of the canonical and Hamilton-Jacobi equations 


With respect to an orthonormal basis of the configuration space C, the momentum 

multivector may be expanded as (see Eq. (|A20)) 


p= E 

\.J\=D 


Pj = P- ej. 


(Bl) 


The surface element dT can be thought of as an outer product of D infinitesimal vectors, 


dr = dqi A ... A dg^,. 


(B2) 


where each dqj can be decomposed as dqk = dqj.ej (summation over the repeated index from 1 to 


D + N is implied). Denoting the derivative in direction Cj by JP = e 


derivative in direction ej = Cj^ A ... A e 
(131 are cast equivalently as follows: 




9Pj — 


dqj 


j ■ dq, and the multivector 


dp, the canonical equations of motion 


^ ■ (dqi A ... A dqp) = det(dgf), 

^ dH{qi,Pi) \dqi-dqP- Cj = dq^^ for D = 1 

dqj I (dr • dq) ■ {Cj ■P) = {§P- Cj) ■ {Ck ■ dr) for D > 1, 

id(g„P/) =0. 


(B3a) 

(B3b) 

(B3c) 
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The determinant in the first equation arises from Formula (A16|. The subsequent Formula (A17) 
can be employed to cast the second equation, case H > 1 , as 


dH{q, Pi) ^ ^ ((ej^ A ... A e^-J • Cj) ■ (e^ • (dqi A ... A dq^)) 

\J\=D 

f)Pj ^ 

= 7r~ XI ^9mminor(dgMZ,m)(-l)'+'"(5yj,. (B4) 

|J|=D l,m=l 


If we parametrize the surface 7 by H coordinates ti, ... ,td, the infinitesimal tangent vectors 
are expressed 


dq 


dqk 


dq^ = Cfc • dqj = ei • v ;—drj = ——drj (j not summed over) 

OTj OTj 


(B5) 


where g = (?(ri,..., td) is the embedding mapping from the parameter sp ace to the configuration 
space. The corresponding Jacobian appears on the right-hand side of Eq. (B3a), which now reads 


det {dqf*‘) = det 


dqj^ 


'3 k 


dn 


dri .. .dr, 


D, 


(B 6 ) 


while the infinitesimal element dri... dru can be divided out to renormalize the Lagrange mul¬ 
tiplier A. 


Let us now focus on the local Hamilton-Jacobi equation (20 1 . In components we have 

H{qi, {dq A S) ■ ei) = 0, 


(B7) 


with 


{dq AS )-61 = (-1)^"^^ • (efe • (e*i A ... A e,^)) = ^(-1)^+^ 

* j=i 


D 


dS 




dq^^ 


(B 8 ) 


where S^^ = S ■ A ... A A ... A are the components of a I? — 1-vector S. For 

example, when D = 2, 


{dq AS) - 61 =^ 

dqi2 


dS,, 

%i ' 


(B9) 
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